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■ ABSTRACT 



. Context. Gamma-ray bursts (GRBs) are attractive objects for constraining the nature of dark energy in a way complementary to other 

O' cosmological probes, especially at high redshifts. However, the apparent magnitude of distant GRBs can be distorted by the gravitational 
lensing from the density fluctuations along the line of sight. 

Aims. We investigate the gravitational lensing effect on the cosmological parameters and dark energy equation of state from GRBs. 
(-H Methods. We first calibrated the GRB luminosity relations without assuming any cosmological models. The luminosity distances of 
O ^ low-redshift GRBs were calibrated with the cosmography method using a latest type la supernova (SNe la) sample. The luminosity distances of 
I high-redshift GRBs were derived by assuming that the luminosity relations do not evolve with redshift. Then we investigated the non-Gaussian 
O nature of the magnification probability distribution functions and the magnification bias of the gravitational lensing. 

^-J Results. We find that the gravitational lensing has non-negligible effects on the determination of cosmological parameters and dark energy. 

^ The gravitational lensing shifts the best-fit constraints on cosmological parameters and dark energy. Because high-redshift GRBs are more 

I I [ likely to be reduced, the most probable value of the observed matter density SIm is slightly lower than its actual value. In the ACDM model, 

we find that the matter density parameter CIm will shift from 0.30 to 0.26 after including the gravitational lensing effect. The gravitational 

lensing also affects the dark energy equation of state by shifting it to a more negative value. We constrain the dark energy equation of state out 

, to redshift z ~ S using GRBs for the first time, and find that the equation of state deviates from ACDM at the Itr confidence level, but agrees 

^—T ' with w = -1 at the 2cr confidence level. 
'NT 

o 



Key words, gamma rays: bursts — cosmology: cosmological parameters 

,—1 , 1. Introduction (Aldering et al. 2004). They cannot provide any information 

L' ■ on the cosmic expansion beyond redshift 1.7. On the other 

■ Our Universe is currently in accelerating expansion. This is , , i. , /r^n-n \ j i j j . r .u- 

. ^ ■' at- hand, gamma-ray bursts (GRBs) are ideal candidates tor this 

revealed by observations of high-redshift type la supernovae . . . t-, , . , , . c r^nn i .u j 

■' jff investigation. The high luminosities ot GRBs make them de- 

i2 (SNe la) (Riessetal. 1998; Perlmutteretal. 1999), cosmic mi- . . ,1 . » .i. j ^ .i. • -ui ■ o r u 

M , , ,^^r, • /o , , tectable out to the edge ot the Visible universe (Bromm & Loeb 

Ci , crowave background (CMB) fluctuations (Spergel et al. 2003, 2002, 2006). The farthest GRB is GRB 090429B at z = 9.4 

2007), and the large-scale structure (LSS) (Tegmai'k et al. (Cucchiai-a et al. 2011). Schaefer (2007) compiled 69 GRBs 

2006). These observations suggest that the composition of the , • ,^ , ^ 1 • -4. ■ j- ^ i_- u 

'''' ^ to simultaneously use five luminosity indicators, which are re- 

universe may consist of an extra component such as dark en- j^^.^^^ _ ^ ^^^^.^^ ^ g^^^^jj 2000)^ 

ergy, or that the equatioris governing gravity may need a vari- ^^^^-^^^^ & Ramirez-Ruiz 2000), £peak - L (Schaefer et al. 

ation to explain the acceleration ot the universe at the present ^rvrvox 7- //-^u- 1 j r^i.- n- • o r ^rvrv/i\ j 

^ ^ 2003), £peak - £y(Ghirlanda, Ghisellini & Lazzati 2004), and 

^ ■ trt - L (Schaefer 2007). Here the time lag (ria„) is the time 

In order to measure the expansion history of our Universe, . f^u » u a a t^t- r ■ 1 1 • 

^ ■' shitt between the hard and sott light curves, L is the peak lumi- 

we need the Hubble diagram of standard candles. Type la ^ r^nn • u i-. f u . j . u 

° -'^ nosity ot a GRB, the variability V ot a burst denotes whether 

supernovae that have played an important role in constrain- ,. , , • 1 j u • j u 

^ 1- r jj-g light curve is spiky or smooth and it can be obtained by 

ing cosmological parameters are the well-known standard can- 1 w i- j c u j i- 

° tor calculating the normalized variance ot an observed light curve 

dies. Unfortunately, it is difficult to observe SNe la at z > j j • e ^ ^■ /c • e 

around a smoothed version ot that light curve (Fenimore & 
1.7, even with excellent space-based projects such as SNAP • r. • onnrw r- • »u u * . u- u 

J Ramirez- Ruiz 2000), tpeak is the photon energy at which the 

Send ojfprint requests to: F. Y. Wang (fayinwang@nju.edu.cn) and Z. spectrum peaks, = (1 - cos ^jO^y^so is the collimation- 

G. Dai (dzg@nju.edu.cn) corrected energy of a GRB, and the minimum rise time (trt) in 
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the gamma-ray light curve is the shortest time over which the 
light curve rises by half of the peak flux of the pulse. Qi & Lu 
(2010) found a new correlation in the X-ray band of GRB after- 
glows. After calibrating with luminosity relations, GRBs may 
be used as standard candles to provide information on the cos- 
mic expansion at high redshift and, at the same time, to tighten 
the constraints on cosmic expansion at low redshift (Dai et al. 
2004; Ghirlanda et al. 2004; Friedman & Bloom 2005; Liang & 
Zhang 2005, 2006; Wang & Dai 2006; Schaefer 2007a; Wright 
2007; Wang, Dai & Zhu 2007; Wang 2008; Qi, Wang & Lu 
2008a,b; Liang et al. 2008; Amati et al. 2008; Cardone et al. 
2009; Liang et al. 2009; Qi, Lu & Wang 2009; Izzo et al. 2009). 
Gamma-ray bursts can also potentially probe the cosmographic 
parameters to distinguish between dark energy and modified 
gravity models (Wang, Dai & Qi 2009a, b; Vitaghano et al. 
2010; Capozziello & Izzo 2008; Xia, et al. 2011). Samushia 
& Ratra (2010) also derived constraints on the ACDM and 
XCDM models using a smaller set of 69 GRBs analyzed in two 
different ways, following Schaefer (2007) and Wang (2008). 
They found that GRB data disfavor standard ACDM at a Icr 
confidence level, and is consistent with this model at a 2cr 
confidence level. They also found that the two techniques give 
somewhat different cosmological constraints. This means that 
model-independent caUbration method and more GRB data are 
needed to improve this situation. We will use the cosmographic 
parameters to calibrate the GRB luminosity correlations. More 
recently, Wang, Qi & Dai (2011) enlarged the GRB sample 
and put constraints on cosmological parameters and equation 
of state of dark energy. In this paper, we use this GRB sample 
which includes 1 16 GRBs. 

However, the gravitational lensing by random fluctuations 
in the intervening matter distribution induces a dispersion in 
GRB brightness (Oguri & Takahashi 2006; Schaefer 2007), 
degrading their value as standard candles as well as SNe la 
(Holz 1998). Gamma-ray bursts can be magnified (or reduced) 
by the gravitational lensing produced by the structure of the 
Universe. The gravitational lensing has sometimes a great im- 
pact on high-redshift GRBs. First, the probabiUty distribution 
functions (PDFs) of gravitational lensing magnification have 
much higher dispersions and are markedly farther from the 
Gaussian distribution more remarkably (Valageas 2000a; Wang 
et al. 2002; Ogri & Takahashi 2006). Second, there is effec- 
tively a threshold for the detection in the burst apparent bright- 
ness. With gravitational lensing, bursts just below this threshold 
might be magnified in brightness and detected, whereas bursts 
just beyond this threshold might be reduced in brightness and 
excluded. Schaefer (2007) considered the gravitational lensing 
biases and the Malmquist biases of 69 GRBs. He found that the 
gravitational lensing and Malmquist biases are much smaller 
than the intrinsic error bars. Our method differs in two ways 
from the one used in Schaefer (2007). First, we use the lat- 
est GRB sample that includes 1 16 GRBs. Second, we calculate 
the distance dispersions from the universal probability distribu- 
tion function (UPDF) of the gravitational lensing amplification 
(Wang et al. 2002; Wang 2005). But Schaefer (2007) used the 
method of Gonzalez & Faber (1997), which relies on the poorly 
known luminosity function and number densities of GRBs. 



We explore the gravitational lensing effects on constraints 
of cosmological parameters and equation of state from GRBs. 
We focus on the non-Gaussian nature of magnification proba- 
bility distribution functions and the magnification bias of the 
gravitational lensing. The structure of this paper is arranged 
as follows: in section 2 we caUbrate the luminosity relations 
of GRBs in a cosmological model-independent way. The con- 
straints on the cosmological parameters and dark energy in- 
cluding gravitational lensing are presented in section 3. In sec- 
tion 4 we present model-independent constraints on the dark 
energy equation of state to z ~ 8 including the weak gravita- 
tional effect. In section 5 we summarize our findings and give 
a brief discussion. 



2. Calibration of the luminosity relations of GRBs 

In this section we calibrate the GRB luminosity relations with 
cosmographic parameters. The cosmographic parameters are 
cosmology-independent. The only assumption is the basic sym- 
metry principles (the cosmological principle) that the universe 
can be described by the Friedmann-Robertson- Walker metric. 

The expansion rate of the Universe can be written in terms 
of the Hubble parameter, H = aja, where a is the scale fac- 
tor and a is its first derivative with respect to time. Because we 
know that q is the deceleration parameter, related to the sec- 
ond derivative of the scale factor, j is the so-called "jerk" or 
statefinder parameter, related to the third derivative of the scale 
factor, and s is the so-called "snap" parameter, which is related 
to the fourth derivative of the scale factor. These quantities are 
defined as 



s = 



1 a 
1 a 
J_a 



(1) 



(2) 



(3) 



The deceleration, jerk and snap parameters are dimensionless, 
and a Taylor expansion of the scale factor around provides 

a{f) = ao |l + Hoit - to) - ^qoH^it - tof + ^jWqC^ - ^o)^ 

and so the luminosity distance 

c f 1 1 

^^L = — |z -H -(1 - qo)z^ - -^{^-qo-3ql + jo)z^ 

+ ^ [2 - 2^0 - 15^0 - 15^0 + 570 + lO^oJo + ^oj + 0(z^)}(5) 

(Visser 2004). Cattoen & Visser (2007) pointed out that it is 

useful to recast di as a function of an improved parameter y - 
z/(l + z) and constrained the cosmographic parameters using 
SNe la data. In this way, because z e (0, oo) mapped into y e 
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(0, 1), the luminosity distance with the fourth order in the y- 
parameter becomes 

ddy) = ^{y - {(qo - 3)/ + i[l2- 5qo - (jo + ^o)] / 

+^ [60 - 7;o - lOQo - 32^0 + lO^ojo + ^qo^o + 21ql-l5ql + so] 

+0(y')}, 

where Qq = 1 + kc^ I H^a^(to) is the total energy density. 

We used the Union2 (AmanuUah et al. 2010) dataset to 
fit the cosmographic parameters. We set the Hubble parame- 
ter Hq = 70 km/s/Mpc. The likelihood function for qo, jo, so 
can be determined from;^'^ statistics 




X' 
where 



(7) 



th\2 



, .data , ,th i 

(T ^ err 



C7 



In the calculation, we used Markov chain Monte Carlo tech- 
niques. A Markov chain with samples on the order of 10^ is 
generated according to the likelihood function and then prop- 
erly burned in and thinned to derive statistics of the parameters 
of interest ^o, Jo and so- The best-fitting results are 

qo = -0.46±0.15 , jo = -1.85±2.25, so = -25.53+35.30.(9) 

The probability distribution of qo, jo and so are shown in Fig.[T] 
Our results are consistent with those in Vitagliano et al. (2009), 
but our results are tighter. Gao et al. (2010) and Capozziello & 
Izzo (2010) used the Union dataset to calibrate the correlations 
of GRBs. Izzo, Luongo & Capozziello (2010) also used the 
Union2 dataset to constrain the cosmography parameters. 





Table 1. Calibration results of five GRBs luminosity relations. 



relations 


a 




b 




L - Tlag 


52.01 ± 


0.19 


-0.83 ± 0.02 


0.52 


L-V 


51.29 ± 


0.32 


0.85 ± 0.32 


0.84 


^ ~ ^peak 


51.79 ± 


0.16 


1.41 ±0.06 


0.62 


L - trt 


52.72 ± 


0.23 


-1.48 ±0.26 


0.58 


~ ^peak 


50.63 ± 


0.10 


1.62 ±0.04 


0.18 



Fig. 1. The probability distribution of qo, jo and so. 



relations, we used the same method as Schaefer (2007) to cal- 
culate the real distance modulus. 



We used the latest GRB sample and five luminosity rela- 
tions in Wang et al. (2011). We first used the best-fitting cos- 
mography parameters (see Eq|9]l to calculate the luminosity 
distance of GRBs at z < 1.4. Then we fitted the five lumi- 
nosity relations at z < 1 .4. The calibrated relations are shown 
in Tabl^U The intrinsic scatter of the L-V relation is too 
large, therefor we discard it below. The calibrated luminosity 
relations are completely cosmology-independent. Wang et al. 
(201 1) established that these relations do not evolve with red- 
shift and are valid at z > 1 .40. The luminosity or energy of a 
GRB can be calculated at high redshifts. In this way the lumi- 
nosity distances and distance modulus can be obtained. After 
obtaining the distance modulus of each burst with one of these 



Hii,^(Y,l^i/(rl)/(Y^CT-^f), (10) 

where the summation runs from 1-4 over the relations with 
available data, is the best-estimated distance modulus from 
the i-th relation, and cr^, is the corresponding uncertainty. The 
uncertainty of the distance modulus for each burst is 

^,rn-i^^,!r'^'- (11) 

The calibrated GRB Hubble diagram is shown in Fig. [2] 
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Fig. 2. Calibrated Hubble diagram of GRBs. The line shows the 
distance-redshift relation in a flat ACDM (Q^ = 0.3). 




Fig. 3. The values of -/f™„ in the flat ACDM (Qm = 0.3). 

3. Magnification probability distribution of 
gravitational lensing 

Owing to the gravitational lensing, the convergence k is given 
by (Bernai-deau et al. 1997; Kaiser 1998) 



3 r^' 



(12) 



Here;^' is the comoving distance (xs corresponds to the redshift 
Zs of the source). 



Ax = 



Dix) = Sinn 



and where Q,t = I-Qm-^^a- The density contrast (5 = (p-p)/p. 
We can see from Eq.(lO) that there exists a minimum value of 
the convergence: 

3 r^' 

Kmin = - - Qm J Wix, Xs)- (13) 

Accordingly the minimum magnification factor is T,„,„ = 1/(1- 
Kmin)^- The plot of /c„„„ is shown in Fig. [3] We can derive the 
magnification PDF for an arbitrary cosmological model by us- 
ing the universal probability distribution function (UPDF) of 
the gravitational lensing amplification (Wang et al. 2002; Wang 
2005). The UPDF can be fitted to the stretched Gaussian (Wang 
et al. 2002), 

w if I j 



g(T]\^,,) = Cnorm exp 



(14) 



where C„onn, ^Ipeak, w, and q depend on and are independent 
of 77. The parameter 77 is defined by (Valageas 2000a, 2000b) 



J] = 



1 - T„ 



= 1 + 



K _ if d;\fvv(y,;\f,,) (p/p) 

<mm\ jQ'dxw(x,Xs) 



(15) 



where t is the magnification factor Note that 77 is the average 
matter density relative to the global mean. The variance of rj is 
given by (Valageas 2000a,2000b; Wang et al. 2002) 

^'^lo '^^iy) ^^^^ 
Here 



luiz) 



n dk 

Jo k 



^Xwix^Xs), 
A\k,z) 



W\Dk0o), 



where A^(k,z) - A-nk^P{k,z), k is the wavenumber, and P(k,z) 
is the matter power spectrum. We used the non-linear power 
spectrum from Peacock & Dodds (1996), which is based on 
N-body simulations. The matter power spectrum is shown in 
FigS W(Dk6o) = 2Ji(Dk0Q)/(Dk6Q) is the window function 
for smoothing angle Oq. Because GRBs are point sources, we 
adopted a sufficiently small smoothing angle 60 = 0.1" (Oguri 
& Takahashi 2006). Here Ji is the Bessel function of the first 
kind. 

Using T = 1 + 2|/c„„„|(7/ - 1) we can obtain 



p{r) 



(17) 



For an arbitrary cosmological model, one can compute from 
Eq.(14), and then the UPDF and p{T) can be computed. In Fig|5] 
we present the magnification probability distribution functions 
P(t) at redshifts z = 1, z = 3 and z = 7 with Q.m = 0.27 and 
Qa = 0.73. The results agree with the /?(t) derived in Wang 
et al. (2002), Oguri & Takahashi (2006) and Holz & Linder 
(2005). The probability distribution functions at high redshifts 
have a higher variance and a lower height of the maximum. 
The peaks reduce to a smaller magnification factor r. From the 
probability distribution functions, we can see that the gravita- 
tional lensing is crucial for high-redshift objects. 
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Fig. 4. Linear and the nonlinear matter power spectrum at 
z = 1 . We used the matter transfer function from Bardeen et 
al. (1986). The parameters are Qm = 0.27, h = 0.7, Q;, = 0.04 
and (Tg - 0.96 (from the WMAP seven-year results). 



0.8 




0.8 1.0 1.2 1.4 



Fig. 5. Magnification probability distribution functions of grav- 
itational lensing at redshifts z = 1, z = 3 and z = 7. 

4. Constraints on cosmological parameters and 
dark energy including magnification bias 

The random magnification of distant sources by the gravita- 
tional lensing induces some bias in the observed sample. We 
calculated the magnification bias as follows. Once the distri- 
bution p{t) of magnification factor t is computed, the mag- 
nification bias is drawn from the distribution pir). The dis- 
persion of distance modulus by the gravitational lensing is 
Aju - -2.5Logr. In Fig|6]we show the dispersions of the GRB 
distance modulus induced by the gravitational lensing. We can 
derive the likelihood function from Bayes' rule (e.g., Gregory 
2005). For each point in a cosmological parameter space, the 
likelihood function is determined by convolving the distribu- 



Fig. 6. Dispersions of GRBs distance modulus induced by 
gravitational lensing. 

tion of T with the intrinsic dispersion distribution (Dodelson & 
ValHnotto 2006). We used h = 0.742 ± 0.036 in our calculation 
(Riess et al. 2009). 

4.1. The ACDM cosmology 

The luminosity distance in a Friedmann-Robertson-Walker 
(FRW) cosmology with mass density Qm and vacuum energy 
density (i.e., the cosmological constant) Oa is (Carroll, Press 
& Turner 1992) 

dL = c(l+z)i/o'|Q*r'^2sinn{|Q,|'/2 

X r dz[{l+z)\l+nMZ)-z(2 + z)nAr^'^]- (18) 
Jo 

We used the GRB sample to constrain the cosmological pa- 
rameters. FigQ shows the Icr to 3cr contour plotting in the 
Q.M - plane. The black line contours from 1 16 GRBs show 
= 0.30+° '^!^ and Qa = 0.84:|:JJ-^° (Icr). The dashed contours 
from 116 GRBs including gravitational lensing magnification 
bias show D.m = 0.26+JJ-,^l^ and Qa = 0.87+J}26 (i^-) prom the 
two contours we can see that the gravitational lensing biases 
the constraints on cosmological parameters. Because the solid 
line in Fig|7]represents a flat universe, our result agrees with a 
flat universe. The contours of GRBs at higher redshifts are al- 
most vertical to the Qm axis because the cosmology is matter- 
dominated at high redshifts. 

4.2. The w{z) = wq model 

We consider an equation of state for dark energy 

w(z) = wo. (19) 

In this dark energy model, the luminosity distance for a flat 
universe is (Riess et al. 2004) 

^/i = c//o'(l+z) rdz[(l+zfnM+{l-^M)(l+zf^'^""'^V''^-(20) 
Jo 
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0.0 0.2 



3 4 5 

redshift z 



Fig. 7. Confidence contours of likelihood from Icr to 3 cr in the 
ACDM model. The solid-line and dotted-line contours from 
1 16 GRBs without and with the magnification bias. 





redshift z 



Fig. 9. Estimates of the uncorrected dark energy EOS parame- 
ters w(z)- Top: uncorrected dark energy parameters versus red- 
shift, in which the vertical errorbars correspond to Icr and 2cr 
confidence levels of w(z) and the horizontal errorbars span the 
corresponding redshift bins. Bottom: Same as top panel but in- 
cluding the weak lensing effect. 



Fig. 8. Same as Fig|7]but in w(z) - wo model. 



Fig[8]shows the constraints on wq versus in this dark energy 
model. The black solid line contours give constraints from 116 
GRBs and we have Qm = 0.29^°ll(lcr) and wo = -1.1. The 
dashed contours give constraints from 116 GRBs including a 
gravitational lensing magnification bias: D.m = 0.26;^!! 25(10-) 

r+0.80 
-"-2.03- 



and wo = -1.05+° '*° 



5. Model-independent constraints on the dark 
energy equation of state 

We first briefly describe a model-independent method to con- 
strain the equation of state (EOS) (for more details, see Qi, 
Wang & Lu 2008a). We adopt the redshift binned parametriza- 
tion for the dark energy EOS as proposed in Huterer & Cooray 
(2005), in which the redshifts are divided into several bins and 
the dark energy EOS is taken to be constant in each redshift 
bin but can vary from bin to bin. For this parametrization, f(z) 
takes the form (Sullivan et al. 2007) 



f(z„-i <z < z„) = (1 +z) 



3(I+R'„) 



(21) 



!=0 
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where w,- is the EOS parameter in the / redshift bin defined 
by an upper boundary at z,, and the zeroth bin is defined as 
Zq = 0. This parametrization scheme assumes less about the 
nature of the dark energy, especially at high redshift, compared 
with other simple parameterizations, because independent pa- 
rameters are introduced in every redshift range and it could, 
in principle, approach any functional form with increasing the 
number of redshift bins (of course, we would need enough ob- 
servational data to constrain all parameters well). For a given 
set of observational data, the parameters w, are usually corre- 
lated with each other, i.e. the covariance matrix 

C = <ww'^> - <w)<w'^), (22) 

is not diagonal. In the above equation, w is a vector with com- 
ponents iv, and the average is calculated by letting w run over 
the Markov chain. A new set of dark energy EOS parameters 
Wi defined by 

w - Tw (23) 

is introduced to diagonalize the covariance matrix. The trans- 
formation of T advocated by Huterer & Cooray (2005) has the 
advantage that the weights (rows of T) are positive almost ev- 
erywhere and localized fairly well in redshift, which facilitates 
an interpretation of the uncorrected EOS parameters vv,. The 
evolution of the dark energy with respect to the redshift can be 
estimated from these decorrelated EOS parameters. The trans- 
formation of T is determined as follows. First, we define the 
Fisher matrix 

F = C"' = O'^AO, (24) 

where O is orthogonal matrix and A is diagonal. Then the trans- 
formation matrix T is given by 

T = 0'^AiO, (25) 
except that the rows of the matrix T are normalized such that 

j 

In addition to the Union2 SNe la sample and 116 GRBs, we 
also used the distance ratio Dv(035)/Dv(0.2) = 1.736 ± 0.065 
from SDSS7 data (Percival et al. 2010) and the shift parameter 
R - 1.725 ± 0.018 from the WMAP seven-year data (Komatsu 
et al. 2010). The four redshift bins are 0-0.2, 0.2-0.5, 0.5 -1.8 
and 1.8 - 8.5. We marginalize over Hubble parameter Hq, as- 
suming a broad uniform prior over the range 50 < Hq < 85 
km s"' Mpc '. We also marginalize over Om assuming the 
quoted prior from Komatsu et al. (2010). Fig. |9] shows the con- 
straints on EOS w(z). This is the first time that the EOS is 
constrained beyond the redshift 1.7. From this figure we can 
conclude that even though the EOS deviates from A at Icr con- 
fidence level, it agrees with w = -1 at a 2cr confidence level. 

6. Discussion and conclusions 

We have presented the gravitational lensing eff'ects on con- 
straints of cosmological parameters and dark energy from 



GRBs. We mainly focussed on the non-Gaussian nature of 
magnification probability distribution functions and the mag- 
nification bias of gravitational lensing. We first used an SNe 
la sample to calibrate the luminosity relations of GRBs. 
Because the luminosity distances of SNe la are completely 
cosmological-model-independent, the GRB luminosity rela- 
tions can be calibrated in a cosmology-model-independent 
way. Then we calculated the PDFs of gravitational lensing. The 
probability distribution functions at high redshifts have higher 
variance and a lower height of the maximum. The peaks reduce 
to a smaller magnification factor t. From the probability dis- 
tribution functions we can see that the gravitational lensing is 
more important for high-redshift objects. Finally we presented 
constraints on cosmological parameters and dark energy. We 
found that the gravitational lensing had non-negligible efifects 
on the determination of cosmological parameters and dark en- 
ergy. The gravitational lensing shifts the best-fit constraints 
of cosmological parameters and dark energy. Because high- 
redshift GRBs are more likely to be reduced, the most probable 
value of the observed matter density Qm is slightly lower than 
its actual value. The gravitational lensing also biases a more 
negative value of the dark energy equation of state. We also 
constrained the dark energy equation of state out to redshift 
z ~ 8 in a model-independent way using GRBs for the first 
time, and found that the equation of state deviates from ACDM 
at the Icr confidence level, but agrees with w = -1 at a 2cr 
confidence level. 

As shown in Samushia & Ratra (2010), the cosmological 
constraints from the two methods of Schaefer (2007) and Wang 
(2008) may be different when using 69 GRBs. Therefore we 
emphasize that we need detailed studies of new correlations 
with a much greater number of GRBs and an examination of 
systematic errors to be able to regard GRBs as more accu- 
rate standardizable candles. Now ongoing missions like Swift, 
Fermi and Suzaku, and the collaboration of many observers on 
ground will promise the progression of GRB cosmology. 
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Appendix A: Calculation of the lensing power 
spectrum 

The linear matter power spectrum A^(A;, z) is parameterized as 

Alik,z) = Ak"'^^T\k)D\z), (A.l) 

where D(z) = g(z)/il + z)g(Q) is the linear growth factor, T{k) 
is the transfer function, A is the normalization factor and 
is the primordial fluctuation spectrum. We use the Harrison- 
Zel'dovich spectrum n, = 1 throughout. For the ACDM model 
(w = -1), the relative growth factor g(z) is well approximated 
by (Carroll et al. 1992) 



gA(z) 

with 



(5/2)Qm(z) 



n'Jj'iz) - Qa(z) + (1 + f2M(z)/2)(l + f2A(z)/70) 



(A.2) 



Q.Mil + z)^ + Q., 



f^A(z) = 



f^M(l +Zy + Q.A 



(A3) 



For the transfer function, we adopt the fitting result of Bardeen 
et al. (1986) for an adiabatic ACDM model 



ln(l +2.34^) 
2.34^ 



[l+3.89^+(16.1#+(5.46^)"+(6.71^)'^]-'^^(A.4) 



where q = k/hF, and h = //o/(100kms"' Mpc"'),and F = 
Qm/i exp[-Qi(l + y[2hlQ.My\ is the shape parameter with 
baryon density Q^,. 

For the non-linear power spectrum we adopt the formula 
given by Peacock & Dodds (1996), 



kL = [l+AldkNL)r''^kML, 

l+Bfix + (Ax)"'^ 



fNlix) = X 



1 + [(Ax)V(z)/(Vxi/2)p 
The parameters in the non-linear function fyi are 



(A.5) 



A 


= 0.428(1 + nJ3y 


-0.947 


B 


= 0.226(1 -H«,/3)" 


-1.778 


a 


= 3.310(1 -i-n,/3)" 


-0.244 


P 


= 0.862(1 -H«,/3)" 


-0.287 


y 


= 11.55(1 -H«,/3)" 


-0.423 



which are fitted to the numerical simulation results. 



